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NATURAL FREQUENCIES OF CONTINUOUS 
FLEXURAL MEMBERS 


A. 8. Veletsos,* J. M. ASCE, and N. M. Newmark,** M. ASCE 


SYNOPSIS 


A numerical method is presented for calculating the natural frequencies 
of undamped flexural vibration of continuous beams on rigid supports and of 
rigid jointed plane frameworks without sidesway. The method is based on 
concepts which are familiar to structural engineers, and it is reduced to a 
routine system of numerical calculations which, when repeated a sufficient 
number of times, will give the natural frequencies of the system considered 
to any desired degree of accuracy. Numerical values for the various physi- 
cal constants which are necessary in the analysis by this method are tabulated 
in the Appendix. The application of the method is illustrated by three numeri- 
cal examples. 


INTRODUCTION 


The problem of calculating natural frequencies of vibration of structures 
is of great practical importance in the analysis and design of structures sub- 
jected to time-dependent or dynamic forces. This problem is particularly 
significant in the study of the behavior of structures under stationary periodic 
forces, such as those which arise from rotating or reciprocating machinery. 
If the frequency of the disturbing force is sufficiently close to one of the natu- 
ral frequencies of the structure, the structure will undergo vibrations which, 
in the absence of appreciable damping effects, may become prohibitively large. 

The purpose of this paper is to describe a method for calculating the natu- 
ral frequencies of flexural vibration of elastic beams continuous over rigid 
supports, and of rigid jointed, plane, elastic frameworks without sidesway. 
The mass of the members composing the structure is assumed to be uniform- 
ly distributed along each member. A system with distributed mass and elas- 
ticity has an infinite number of natural frequencies. With the method pre- 
sented herein one is capable of determining as many of the natural frequencies 
as one desires. 

The application of classical analytical methods for determining natural 
frequencies of continuous beams and frames has been described, among 
others, by Darnley,(1)1 Timoshenko,(2) Bennon,(3) and Hohenemser and 
Prager.‘ ) Additional methods having a fairly general applicability have been 
presented by Mudrak(5,6) Saibel(7) and Looney \8),A.M. ASCE, and a numeri- 
cal method applicable primarily to fundamental frequencies has been de- 
veloped by Masur. (9 


* Research Asst. Prof. of Civ. Eng., Univ. of Illinois, Urbana, Ill. 
** Research Prof. of Structural Eng., Univ. of Illinois, Urbana, Ill. 
1. Numbers in parentheses, unless otherwise identified, refer to the Bibli- 
ography at the end of the paper. 
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The method presented in this paper is an extension of the well-known 
Holzer’s method(19) for calculating natural frequencies of torsional vibration 
of shafts. It is based on principles which are familiar to structural engineers, 
and it is reduced to a straight-forward system of numerical calculations 
which, when repeated a sufficient number of times, will give any of the natural 
f-equencies of the system considered to any desired degree of accuracy. This 
procedure is simple, of wide applicability, and appears to be less subject to 
computational errors than any of the other methods known to the authors. 

This paper is divided into three parts and an Appendix. Part I describes 
the principles underlying the method and gives definitions of the various 
physical constants entering in the analysis. Numerical values of these con- 
stants are tabulated in the Appendix. The application of the method to con- 
tinuous beams is described in Part II, and the extension of the method to con- 
tinuous frames is explained in Part III]. For convenient reference, a detailed 
outline is included of the steps involved in the application of the method to 
particular problems, and the method is illustrated by three numerical exam- 
ples. 


Part I. General Description of Method 


Basic Principles 


The method is based on the fact that, in the absence of damping, the excit- 
ing couple which is necessary to maintain a dynamical system in a steady- 
state of forced vibration with finite amplitudes becomes equal to zero when 
the frequency of the couple is equal to one of the natural frequencies of the 
system. It is assumed that the exciting couple is applied at a joint for which 
the amplitude of rotation is different from zero for all the natural frequencies 
that are to be determined. 

Consider that a continuous beam or frame is in a steady-state of forced 
vibration under the action of a harmonically varying exciting couple applied 
at a suitable joint, designated as joint z. The amplitude of this couple is as- 
sumed to be such that the amplitude of the rotation (or of the internal bending 
moment) at a joint of the structure different from z, say at joint 1, has a pre- 
scribed finite value. The vibration of the structure is then harmonic and its 
frequency is equal to that of the exciting couple. Furthermore, since the ef- 
fect of damping is neglected, the amplitude of vibration is constant and the 
response is either in phase with or 180° out of phase with the exciting couple. 

For a given system, the magnitude of the exciting couple necessary to 
maintain the prescribed amplitude of vibration at joint 1 depends on the fre- 
quency of the couple. For the limiting case of zero frequency, the magnitude 
of the couple is obviously finite; if desired, its value may be calculated by any 
of the available methods of indeterminate stress analysis. As the frequency 
of the exciting couple increases above zero, the structure is acted upon by 
distributed inertia forces of increasingly greater magnitudes. These forces 
are at first in phase with the exciting couple and produce displacements in 
excess of those produced by the couple acting statically. With a further in- 
crease in the frequency of the exciting couple, the inertia effects predominate, 
and for a frequency equal to the fundamental or lowest natural frequency of 
the structure, the magnitude of the couple becomes equal to zero. This con- 
dition of the vanishing exciting couple is characteristic of all the natura! fre- 
quencies of an undamped elastic system. 

Therefore, to determine the natural frequencies of a structure the follow- 
ing approach may be used: (a) consider that the structure oscillates in a 
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steady-state at an assumed frequency of vibration and with a fixed amplitude 
of rotation (or internal bending moment) at joint 1; (b) determine the magni- 
tude of the exciting couple at joint z necessary to maintain the vibration; 

(c) repeat these steps for a number of assumed frequencies; and (d) plot the 
magnitude of the exciting couple as a function of the frequency of vibration. 
The frequencies for which the exciting couple vanishes represent the desired 
natural frequencies of the structure. 

It should be noted that the amplitude of rotation (or internal bending mo- 
ment) at joint 1 is assumed to be different from zero in this method. Accord- 
ingly, if there exists a natural frequency for which this assumption is not 
satisfied, the method will fail to detect it. Such a condition may occur in the 
case of continuous frames. A supplementary technique for determining the 
natural frequencies which the main procedure may fail to detect is described 
in Part III of the paper. : 

The primary computational task in the method outlined is the evaluation of 
the exciting couple. This may be accomplished by a number of different pro- 
cedures. The conditions to be satisfied are simply those of equilibrium and 
continuity for each joint of the structure. These may be expressed by a set 
of linear equations in a number of different ways, and the equations may be 
solved by a variety of procedures. In the procedure presented herein, these 
conditions are expressed in the form of a dynamic three-slope equation, and 
the distortions of the joints are computed by the repeated application of this 
equation, working progressively from one end of the structure to the other. 


Elastic Constants for a Vibrating Bar 


The following terms refer to a bar which is simply supported at one end 
and fixed at the other, and is in a steady state of oscillation under the action 
of a harmonically varying bending moment applied at the simply supported 
end. 

The amplitude of the bending moment necessary to produce, at the end at 
which the exciting moment is applied, a steady-state rotation of unit amplitude 
is defined as the “dynamic flexural stiffness” and will be denoted by the sym- 
bol K. 

The ratio of the bending moment developed at the far fixed end to the excit- 
ing moment at the near end is defined as the “dynamic flexural carry-over 
factor,” and will be denoted by k. 

The foregoing terms are generalizations of those introduced by Hardy 
Cross,(11) Hon. M. ASCE, for the analysis of frames under static loads, and 
they were first used by R. E. Gaskell (12) in applying the method of moment 
distribution to the analysis of the steady-state forced vibration of frames. 

For certain conditions of symmetry, antisymmetry, and for those cases 
for which the degree of restraint at the far end of a member is known, it 
may be convenient to use effective stiffnesses. The effective stiffness for a 
symmetrically deformed bar is 


K = K(I-k). (1) 
For an antisymmetrically deformed bar, the corresponding stiffness is 


K° = K(I+k). (2) 
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For a bar, the far end of which is subjected to a rotational elastic restraint 
of stiffness R, the effective stiffness is 


2 
(3) 


These expressions are identical to those for the static case. (13,14) 
The discussion in this paper will be restricted to members having uniform 

cross section and mass per unit of length. The stiffness and the carry-over 

factor for both ends of such a member are equal. 


Numerical Values of Constants 


The peitinent analytical expressions for the quantities defined in the pre- 
ceding section can be found in reference (12). The dynamic stiffness may be 
expressed as 


K= (4) 


where Cy is a dimensionless factor depending on the dimensionless parame- 
ter 


(5) 


In the foregoing expressions, m is the mass per unit of length of the bar; w 
is the circular frequency of vibration; E is the modulus of elasticity of the 
material in the bar; I is the moment of inertia of the cross section of the bar 
about its centroidal axis; and L is the span length of the bar. The carry-over 
factor k is dimensionless and depends also on the parameter A. 

A graphica) representation of the variation of the stiffness and the carry- 
over factor as a function of A is given in Fig. 1. For A = 0 these quantities 
have the well known static values of K = 4EI/L and k = 0.5. Values of A equal 
to 3.927, 7.069 and 10.210 correspond, respectively, to the first, the second, 
and the third natural frequencies of a bar fixed at one end and simply support- 
ed at the other. For these frequencies, the value of the dynamic stiffness is 
equal to zero since no exciting moment is required to maintain the vibration. 
Furthermore, since the amplitude of the moment at the fixed end of the bar is 
different from zero for these frequencies, the value of the carry-over factor 
is infinite. Values of A equal to 4.730 and 7.853 correspond to the first and 
the second natural frequencies of a bar fixed at both ends. For these fre- 
quencies the stiffness of the bar becomes infinite since the end rotations are 
zero. 

In the procedure presented herein, only the stiffness and the product of the 
stiffness and the carry-over factor are needed. Numerical values of these 
quantities have been computed to seven significant figures(15) for values of A 
between zero and 10.20 at increments of 0.01. A portion of these results are 
tabulated to five significant figures in the Appendix. It should be noted that 
some of these values could have been obtained from functions tabulated in 
reference (4). 
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Part Il. Continuous Beams 


Characteristics of Beams Considered 


The beams are assumed to be initially straight and continuous over non- 
deflecting supports. The cross section and the mass per unit of length of each 
beam may vary from one span to the next, but in any one span these quantities 
are considered constant. It is assumed that the vibration is restricted to one 
of the principal planes of flexure of the beam, and that the cross sectional di- 
mensions of each span are small in comparison to its length, so that the ef- 
fects of shearing deformation and rotatory inertia are negligible. 


Sign Convention and Nomenclature 


Downward deflections and clockwise rotations are taken as positive. In- 
ternal bending moments in a member at its connection with other members, 
and external couples are taken as positive when acting in a clockwise direc- 
tion, 

The supports are numbered consecutively starting with 1 at the extreme 
left end and terminating with z at the extreme right end. The quantities Lj, 
EjIj, mj, Aj, K,, and (KK); refer to the span between supports j andj +1. The 
symbol 6; denotes the amplitude of the rotation at the j-th support and Mj de- 
notes the corresponding amplitude of internal bending moment. The sub- 
scripts L and R designate, respectively, sections just to the left and just to 
the right of the support. The amplitude of the exciting couple at support z is 


denoted by M,. 


Fundamental Equations 


Figure 2 shows, in an exaggerated scale, the extreme deflected position of 
spans j-1 and j of a continuous beam undergoing steady-state forced vibration. 
The vibration is assumed to be maintained by a single exciting couple applied 
at support z. The rotations and bending moments at the ends of each span are 
indicated in their positive directions. Their value at a time t for support j 


are 


6, (t) = 8;cos wt and M(t) M; cos wt. (6) 


In the equations to be used, the term cos wt appears as a common factor; for 
convenience this will be omitted, and in the remainder of this discussion, the 
terms “amplitude of rotation” and “rotation” and the terms “amplitude of mo- 
ment” and “moment” will be used interchangeably. 

Consider the interior support j. From the conditions of continuity and 


equilibrium it follows that 


(6). (O)r 8; (7) 


(8) 


(M), + = O. 


The moments (Mj);, and (Mj)p can now be expressed in terms of the end 
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FIG. ILLUSTRATIVE EXAMPLE 1 


G 
2 
nso. 
1.35€1 1.35 €1 
™ 0.60™ 1.20™m m 
2 3 4 5 
(a) 
2.40 2.84 2.33 3.34 
0.50 3.665 3.302 3704 2.481 
(1.00) (0.80) @.35) (0.90) 
(b) 
2.893 
0.00) @I.48) 
|_| 
7135-7 


rotations of spans j-1 and j, respectively. To start with, consider span j. 
Assume, first, that the right end of the span is held fixed while the left end is 
rotated through an angle 63; then, the moment at the end being rotated is 
equal to the product of the rotation 6; and the flexural stiffness of the mem- 
ber Kj. Imagine, next, that = left end of the span is kept fixed while the 
right end is rotated through @ +1; the moment induced at the fixed left end is 
equal to the rotation 6; +] multiplied by the product of the flexural stiffness 
and flexural caameanll factor of the member (kK) ;. Since the governing dif- 
ferential equation is linear and the principle of on position is valid, 


(Mp = + (90) 


By considering span j-1, one obtains in a similar manner the result 


(M), = K,.,0, + (kK), (9b) 


Substitution of Eqs. (9a) and (9b) in Eq. (8) yields 
(kK),.,6,., + [K)., + + (kK),@,, = 0. (10) 


It is convenient to rewrite this equation in the form 


8. 


(10b) 


where Kj , the total flexural stiffness at support j, is 


Ky = + (MN) 


The foregoing expressions are, strictly speaking, applicable to interior 
supports only. For the end supports, the equations have the following special- 
ized forms. Consider the general case in which the extreme ends of the beam 
are elastically restrained against rotation. The relationship between the re- 

straining end moments and the corresponding rotations are 


M, = -R,@, and M, = -R,@,, (12) 


where Rj and R, are known stiffnesses of the rotational restraints at the left 
and the right ends, respectively. For a hinged end R = 0, and for a fixed end 
R = infinity. The negative signs in the foregoing expressions denote that for 
a positive restraint (positive value of R), the restraining moment acts in a 
direction opposite to the direction of rotation of the span. 

Consider first the extreme left end of the beam. For a hinged or partially 


fixed condition, Eq. (10a) reduces to 
+ (kK),@, = 0, (13a) 
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where 
K, = R, + K,. (I3b) 


For a fixed condition, 6; = 0, and the relationship between the internal mo- 
ment at the fixed end and the rotation at the second support is 


M, = (kK), @,. (14) 


Consider finally the right end of the beam. For a hinged or partially fixed 
end, the magnitude of the exciting couple is obtained from a specialized form 
of Eq. (10a) as 


M, = (kK),.,9,., + K,®, (15a) 


where 


K, = K,., + R, (15b) 
At a natural frequency M, must be identically equal to zero. For a fixed end 
0, (16) 


and the criterion for a natural frequency is that Eq. (16) be satisfied. The 
magnitude of the internal bending moment at the fixed end is of no practical 
interest but, should it be desired, it can be calculated from Eq. (15) by taking 
62 = 0. 


Outline of the Procedure 


As previously remarked, the principal computational task of the procedure 
is the determination of the exciting force at the right end of the beam neces- 
sary to maintain a given amplitude of rotation (or internal bending moment) 
at the left end. Since the natural frequencies of a system depend on the rela- 
tive rather than the absolute values of the displacements, any arbitrary am- 
plitude, consistent with the actual boundary conditions may be chosen for the 
amplitude of vibration at support 1. For a hinged or a partially fixed end, 6; 
is taken, for convenience, equal to unity; for a fixed end, 6, is equal to zero, 
and Mj or Mj times the L/EI of some reference span r is taken equal to unity 
instead, 

The procedure may now be outlined as follows: 


1, Start by assuming a value of A for some reference span r; this step is 
equivalent to assuming a frequency of vibration 


2 
Ay / 
Lev 
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Compute then the values of A for the remaining spans from the expression 


im, 


The pertinent calculations are conveniently carried out in a tabular form, as 
shown in the following illustrative example. Equation (17) is obtained readily 
by equating the expression for w of the r-th span to that of the j-th span. 


2. From the table in the Appendix, obtain the values of K and kK for the 
various spans. 


3. Compute the rotation of the beam over the second support from Eq. (13) 
or (14) and, by successive applications of Eq. (10b), evaluate the rotations 63 
to 6,. A convenient scheme for arranging the computations is described in 
the illustrative example. If support z is fixed, the determination of the rota- 
tion 6z completes one cycle of the procedure (see Eq. 16). However, if this 
support is hinged or only partially fixed, it is necessary to carry out the ad- 
ditional step of evaluating Eq. (15). 


4. Repeat steps 1 through 3 for different values of A,, and plot the values 
calculated for 67 or M, as a function of A,. The zero intercepts of the re- 
sulting curve which, in general, has the shape shown in Fig. 4, corresponds 
to the natural frequencies of the system considered. 


Illustrative Example 


It is desired to calculate the first five natural frequencies of the four-span 
continuous beam shown in Fig. 3a. The beam is elastically restrained against 
rotation at the left end and simply supported at the right end. The stiffness of 
the end restraint and the properties of the various spans are shown in Fig. 3a. 

For convenience in carrying out the computations, the natural frequencies 
of the system will be expressed in terms of the pertinent properties of a ref- 
erence span r. In this case, we take r = 1. In terms of the EI/L of the r-th 
span, the values of K and kK for a span j are equal to those obtained from the 
Appendix multiplied t by the dimensionless factor 


Eh be 


(18) 


The quantities A j/Ar and aj are evaluated in Table 1. It should be noted 
that these ealoulations are independent of the frequency of vibration. A com- 
plete cycle of the procedure will now be presented for a trial value of A; 

A = 2.40; this value corresponds to a circular frequency of vibration 


(2.40)* 


The values of A for the various spans are recorded at the top of the dia- 
gram in Fig. 3b. These values are obtained as the product of the assumed 
Ar and each of the values in column 5 of Table 1. The values of @ for the 
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various members are recorded in parentheses on the same diagram. The 
numbers above the values of a give the stiffnesses of the various spans, 
whereas those below the values of @ give the corresponding products of the 
stiffnesses and the carry-over factors. These quantities are expressed in 
terms of the EI/L of the span to which they refer, and they are obtained di- 
rectly from the Appendix. 

Before using them in the fundamental e uations, it is convenient to express 
the quantities K and kK in terms of a common factor. On the diagram in 
Fig. 3c, these ‘quantities are given in terms of the EI/L of the reference span. 
On this diagram, the number above each support represents the value of K 
for that support, and the number at the middle of each span represents the 
value of kK for that span. The value of K for the j-th support is determined 
by taking the sum of the products of the ao values and the K values in Fig. 3b 
for spans j-1 and j. Thus, the total stiffness at support 3 is 


K, = + | FE 7.642 
Similarly, the values of kK are obtained by multiplying the values of kK given 
in Fig. 3b by the corresponding values of a. 

With the foregoing information established, one can compute readily the 
rotations of the beam. On the diagram in Fig. 3c, the number in parentheses 
below each support represents the value of the rotation at that support. The 
rotation 61 is taken equal to unity. The rotation 69 is evaluated from Eq. 
(13a) which, it should be recalled, is a specialized form of Eq. (10b). Thus, 


4.165(1.00) _ _ 
0, = -00) 1.846. 


The rotations 63, 64 and 65 are determined by the application of Eq. (10b) to 
joints 2, 3 and 4, respectively. Thus, 


2.256(1.00) + 6.307(-1.846) _ , 
=- 5-035 4.617, 
_ 2. 033(-1.846) + 7.642(4. 617) _ 
6, = 


+ 
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3.004(4.617) + 7-233(-10.50) _ 


2-095 


These results are recorded on the diagram as they are computed. 

The exciting couple, evaluated at the bottom of the diagram in Fig. 3c by 
use of Eqs. (15), is found to be equal to 17.54 EI/L. 

Since Mg is different from zero, the assumed value of A, = 2.40 does not 
correspond to a natural frequency of the system. The physical significance 
of the computed value of 17.54 EI/L is as follows: the negative of this value 
divided by the rotation 65, 


17.54 EI _ EI 


represents the stiffness of a rotational constraint which, if applied at the 
right end of the beam, would make the assumed frequency correspond to a 
natural frequency of the system. 

By repeating several such cycles of computation for different values of 
A1, the curve shown in Fig. 4 was obtained. The zero intercepts of this curve 
correspond to the desired natural frequencies. The first five critical values 
of A are recorded in the figure. The corresponding circular natural frequen- 
cies Wy, in radians per second, are 


6.27 [er [EE 
16.9 EI 

, (a), = =, /= 


The natural frequencies, in cycles per second, may be computed from the 
relation 


= (19) 


In this example, and in the examples to follow, the calculations have been 
carried out to four significant figures. In the particular cycle described, this 
accuracy was greater than actually required. In general, for the determina- 
tion of the lower natural frequencies of continuous beams slide-rule accuracy 
will prove satisfactory. However, for the higher natural frequencies, particu- 
larly in the case of frames involving closed panels, the computations may in- 
volve small differences between large quantities, and it may be necessary to 
retain a minimum of four significant figures. 


Part III. Continuous Frames Without Sidesway 


The following discussion applies to rigid jointed, plane, elastic frameworks 
for which the joints do not translate. The individual members of each frame 
are considered to be of uniform cross section and mass per unit of length. 


735-12 


2k .0 EI 

(ay), = 


MOMENT Ms IN TERMS oF &L 


FIG. 4 VARIATION OF EXCITING MOMENT AS 
A FUNCTION OF A;, EXAMPLE 1 
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The vibrations are assumed to take place in the plane of the framework. The 
change in length of the members due to axial deformation and the effect of the 
axial forces on the bending moment in the members are neglected. In addi- 
tion, no account is taken of the influence of axial vibrations. As before, the 
cross sectional dimensions of the members are considered to be small in 
comparison to their length, so that the effects of shearing deformation and 
rotatory inertia are negligible. 

For the purpose of this discussion structural frames will be divided into: 
(a) open frames, including those which do not involve any closed panels and 
have known conditions of restraint at all terminals, and (b) closed frames, 
including those which do involve closed panels. Each group will be discussed 
separately. It is first necessary, however, to develop the appropriate slope 
equations applicable to frames. 


Fundamental Equations 


Figure 5 shows s members of a structure rigidly connected at their com- 
mon intersection 0. The ends of the members are assumed to be fixed 
against translation, but free to rotate subject to the restraint imposed by the 
adjoining members. Assume that the structure is in a steady-state of forced 
vibration under the action of an exciting moment applied at a joint different 
from joint 0. 

Let 60,j denote the amplitude of rotation at end 0 of member (0,j), and 
9; g denote the amplitude of rotation at end j of the same member. Similar- 
ly, let Mo Jj and M;, go be the corresponding moment amplitudes at the same 
ends. Also, let Ko. J and (kK) 0 J be, respectively, the v-lues of K and kK for 
member (0,j). 

Since all members are rigidly connected at their joints 


s ® 2 ® “eee 8, = “ef bad (200) 


* 9. (20b) 


Furthermore, since no external moment acts at joint 0, 


Mo, + Mo2 Mo; Mos s O. 
(21) 


The moment Mo ,; may now be expressed in terms of the end rotations of 
member (0,j) as follows: 


Mo) * Koi + (kK), (22) 


Substituting this expression into Eq. (21), one obtains Eq. (23a) 
+ @; = 0, (23a) 
j* 
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which expresses the conditions of both equilibrium and continuity for joint 0 
of the structure. If only two members meet at joint 0, Eq. (23a) reduces to 
Eq. (10a) for continuous beams. 

If the restraints at the far ends of the members meeting at a joint are 
known, it is convenient to use effective stiffnesses. Assume that the re- 

straints at ends 1 and 2 of the portion of the structure shown in Fig. 5 are 

known. Let KO, 1 and Kj}, g be the effective stiffnesses of members (0,1) and 
(0,2), respectively; these may be determined from Eq. (3). Then, Eq. (23a) 
may be written as 


(Ko, + + + 0 
(23b) 


Open Frames 


Simple L-frames and portal frames, such as those shown in Figs. 6a and 
6b, behave as continuous beams on rigid supports. Therefore, their natural 
frequencies can be calculated by the procedure presented in Part II of this 
paper. 

The natural frequencies of open frames, such as those shown in Figs. 6c 
through 6e, can be evaluated by a variation of the procedure applicable to 
beams. Consider any of the frames in Figs. 6c through 6e. Let 1 denote the 
joint of the frame at the extreme left terminal and z denote the joint at the 
extreme right terminal. Without loss of generality, it may be assumed that 
joint z is either hinged or elastically restrained. A fixed end may be treated 
in the same manner as a fixed end of a continuous beam. Assume that the 
structure is in a steady-state of forced vibration under the action of an ex- 
citing couple M, applied at joint z, The amplitude of the slope or of the bend- 
ing moment at joint 1 is assumed to have a fixed value, If the joint is hinged 
or is elastically restrained, the amplitude of the slope is taken equal to unity. 
If the joint is fixed, the amplitude of the internal bending moment is taken 
equal to unity instead. For an assumed frequency of vibration, it is general- 
ly possible to calculate the magnitude of the exciting moment M, in a manner 
entirely analogous to that used for continuous beams, by working progressive- 
ly from one end of the frame to the other, with the only difference that one 
must use the more general Eqs. (23) rather than Eqs. (10). By repeating this 
procedure for several frequencies of vibration, the magnitude of M, may be 
plotted as a function of the frequency of vibration. The frequencies for which 
the amplitude of the exciting couple vanishes are the natural frequencies of 
the frame. 


Natural Modes of Partial Vibration 

When a continuous beam is in a natural mode of svat, each span oscil- 
lates with finite amplitudes. In the case of frames, however, it is possible 
for only a portion of the frame to oscillate while the rest remains stationary 
(vibrates with zero amplitude), Figure 7 presents several such modes for 
the structures shown in Figs. 6c through 6e. Such modes can, of course, 
exist only if the dimensions of the various members satisfy certain definite 
relations. Natural modes for which only a portion of a structure vibrates 
while the rest remains still will be referred to as “natural modes of partial 
vibration.” 

Because of the assumptions on which it is based, the procedure described 
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thus far may fail to detect natural frequencies corresponding to natural 
modes of partial vibration. It may be recalled that the amplitude of the slope 
or bending moment at joint 1 of the structure is assumed to be fixed in the 
procedure. Accordingly, the procedure will fail to reveal those natural fre- 
quencies for which bar (1,2) remains stationary. A second assumption im- 
plicit in the procedure is that the rotation of joint z is different from zero. 
Since this condition is not always satisfied, the procedure will fail also to re- 
veal the natural frequencies for which the bars meeting at joint z are station- 
ary. 

The natural frequencies which the main procedure fails to detect may be 
determined by (a) calculating the frequencies for which bar (1,2) and the bars 
meeting at joint z can remain still, and (b) ascertaining whether or not these 
frequencies are natural frequencies of the system. The details of this sup- 
plementary procedure are described by a numerical example in the following 
section. 

It is worthy of note that natural modes of partial vibration occur rather 
infrequently; they are most likely to occur in the case of frames involving 
panels of uniform members, It is worthy of note also that for most practical 
applications, it may be entirely unnecessary to calculate natural frequencies 
corresponding to such modes. In practice, one is usually interested in deter- 
mining natural frequencies comprised within specified ranges of frequencies. 
Therefore, if the values calculated in the first step of the supplementary pro- 
cedure are found to lie outside the ranges of interest, it will be superfluous 
to carry out the second step. The first step of this procedure, as will be seen 
shortly, can usually be carried out almost by inspection. It should finally be 
noted that natural modes of partial vibration correspond always to higher 
natural frequencies. Therefore, no consideration need be given to these 
modes, if only the fundamental natural frequency of a frame is to be deter- 
mined, 


Illustrative Example 
As an illustration, we shall calculate the first few natural frequencies of 
the simple frame shown in Fig. 8a. The bars identified by (1) are taken 
identical, while the bar designated by (2) is considered to have such dimen- 
sions that 


L, 0.80 


3.927 


The subscripts 1 and 2 refer to bars (1) and (2), respectively. These propor- 
tions have been chosen so that the frame can vibrate in a natural mode of par- 
tial vibration. 

For the sake of brevity, only one cycle of the procedure is presented, The 


computations are given for a value of A; = 3.30, which corresponds to a value 
of Ag = 2.74, 


The appropriate values of K and kK, obtained from the Appendix, are: 


for bars (1), K, = 2-572 E,1,/L, and (kK), = 3.136 g,1,/L ; 
for bar (2), K, = 3-405 g,1,/L, and («K),, = 2.459 
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The data necessary in the analysis are compiled in Fig. 8b. The number 
without parentheses opposite each joint gives the sum of the stiffnesses of 
the members meeting at that joint. The number at the middle of each mem- 
ber gives the product of the stiffness and the carry-over factor for the mem- 
ber. Both quantities are expressed in terms of E41;/ Lj. The numbers in 
parentheses denote the rotations of the joints. These quantities are evaluated 
in the manner described in the following, and are recorded on the diagram as 
they are computed. 

The procedure is started by taking M, = 1.00 E41/ Ly. Then 62 is com- 
puted by application to joint 1 of Eq. (22), 


1.00 _ 
6, = 0.3187 . 


The rotations 64 and 65 are then determined successively by application of 
Eq. (23a) to joints 2 and 4, as follows 


3.138 (0) (0.3187) =-0.522h, 
6, 3138 (0.3187) + 218 + 8.549 (-0.5224) _ 1.409. 


The magnitude of the exciting moment is determined from Eq. (22) as 


[> (1-409) + 2.459 ( 4) | 
« + 2.459 (-0.522 = 5.515 ° 
5 


Since Mg is different from zero, the assumed frequency does not corre- 
spond to a natural frequency. By repeating this procedure for several values 
of Ay, the curve shown in Fig. 9 was obtained, The values of Ay correspond- 
ing to natural frequencies are recorded on the figure. 

It is next necessary to investigate whether there are any natural frequen- 
cies for which My = 0. To this end, the following reasoning is used. In order 
for Mj, to be equal to zero, bar (1,2) must be stationary; then, the rotation 
and the internal bending moment at joint 2 must be also equal to zero. This 
condition requires that bar (2,4) be stationary and that 64 = M42 =0. But 
with joint 4 remaining fixed against rotation, bar (3,4) can oscillate freely 
only at frequencies represented by values of 


= 4-750, 7-853, (24) 


During a natural frequency, every member of a structure vibrates with the 
same frequency. Therefore, in order for these frequencies to be natural 
frequencies of the frame, they must be also natural frequencies of the re- 
maining portion of the frame (bar (4,5) in this particular case). The natural 
frequencies of bar (4,5), considering that its left end is fixed, correspond to 
values of 


ho = 3-927, 7.069, ... (25a) 


which, in turn, correspond to values of 
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A, = 4.730, 8.514, (25d) 


Comparing the values given in Eq. (25b) with those given in Eq. (24) one 
concludes that, within the range of frequencies considered in Fig. 9, the 
values A; = 7.853 and A, = 8.514 do not correspond to natural frequencies, 
whereas the value A, = 4.730 does, The natural mode of partial vibration for 
A 1 = 4.730 is shown by the solid curve in Fig. 10a. 

Finally, it is necessary to investigate whether there are any natural fre- 
quencies for which @5 = 0. One reasons that 65 can be equal to zero only if 
bar (4,5) is stationary, Under this condition, 64 = Mgq 5 = 0; then, bar (3,4) 
can vibrate freely only at frequencies represented by the values of A1 given 
in Eq. (24). Since members (1,2) and (2,4) of the remaining portion of the 
frame are identical to member (3,4), each of the members composing the 
panel (1,2,3,4) can vibrate in a fixed-ended condition for the same frequency; 
therefore, the Ay values given in Eq. (24) correspond to natural frequencies 
of the frame. The natural mode of partial vibration corresponding to Aq = 
4.730 is shown by the dotted curve in Fig. 10a, and the mode corresponding 
to Ay = 7.853 is shown in Fig. 10c. 

It should be noted that a total of three natural modes have been determined 
for a frequency corresponding to a value of A, = 4.730. It can be shown, how- 
ever, that of these three modes only two are independent, the mode shown in 
Fig. 10b being a linear combination of the other two, 

From the foregoing, one concludes that, within the range of frequencies 
considered in Fig. 9, the complete set of Ay values corresponding to natural 
frequencies is 


3.59, 4.22, 4.73 (double),6.80, 7.44, 7.85, and 8.35 


Closed Frames 


The fundamental features of the procedure for calculating natural frequen- 
cies of closed frames are described in this section by reference to the two- 
celled rectangular frame shown in Fig. lla, 

The procedure is started by assuming that the frame is cut at a convenient 
joint. For the example considered, the cut is introduced at joint z = 6, as 
shown in the figure. It is next assumed that the frame undergoes a steady- 
state of forced vibration with finite amplitudes and known frequency under the 
action of an exciting couple applied at the cut joint (z = 6). In general, there 
will be a discontinuity in slope in the cut joint. However, for a natural fre- 
quency of vibration, both the amplitude of the exciting moment and the dis- 
continuity in slope will vanish. Therefore, the criterion for a natural fre- 
quency is that equations 


= (260) 


be satisfied simultaneously. It is, of course, assumed that the frame does 
not oscillate in a natural mode of partial vibration. 


Determination of the Discontinuity in Slope and of the Exciting Moment 

The discontinuity in slope and the magnitude of the exciting moment may 
be determined by working progressively from joint to joint across the struc- 
ture, in the same general way as for open frames, except for one major 
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difference. Whereas for the cases considered thus far the equation for the 
first joint involved only two unknowns, in the present case, the corresponding 
equation involves a minimum of three unknowns. 

To illustrate this point, reference is again made to the frame shown in 
Fig. lla. One starts by applying Eq. (23a) to joint 1. Then, the resulting 
expression will involve the three unknowns, 61, 62, and 93, only one of which 
can be taken equal to a constant. Thus, it iealnes necessary to express 63 
in terms of both 6, and 69. Next, one selects consecutive joints in such an 
order that each new equation involves only one additional unknown, which is 
also expressed in terms of both 6; and @9. The rotation 64 is determined 
by applying Eq. (23a) to joint 2, The rotations 05, 66 4 and 96, 5 are deter- 
mined by the successive application of the same equation to joints 3, 4 and 5, 
respectively. Finally, the exciting moment Mz (z = 6) is obtained from Eq. 
(22). Having expressed all these quantities in terms of 6; and 99, one may 
write Eq. (26a) in the form 


+ = O, 
5,8, + O, (260) 


In these expressions, the c’s are numerical constants, the magnitudes of 
which depend on the assumed frequency of vibration and on the properties of 
all the members of the structure. Equations (26b) are a set of linear, homo- 
geneous equations for the unknowns 6 and 62, and have a solution different 
from zero when the determinant of the coefficients of 9; and 6» is zero: 


4,(6,8,) = (27) 


Co, 


It should be noted that, in this procedure, the cut is always introduced at 
the joint, the rotation of which is evaluated last by the scheme described, 
Furthermore, the number of members meeting at the cut are always equal 
to the number of unknowns necessary to express the rotations of the joints. 

The foregoing criterion of the vanishing determinant will fail to reveal the 

natural frequencies for which 6; and 65 are simultaneously equal to zero. 
It will fail to reveal also the sd Toor traquencies for which the bars meeting 
at the cut are stationary. If necessary, the natural frequencies correspond- 
ing to these conditions of partial vibration may be calculated by the supple- 
mentary procedure described previously. 

Consider now the frames shown in Figs. 11b and llc. The natural fre- 
quencies of vertical vibration of the single-bay, multi-story building frame 
shown in Fig. 11b may be found in identically the same manner as that de- 
scribed for the frame in Fig. lla. The rotations may be expressed in terms 
of 6, and 65, and the cut may be introduced at joint 9 or 10. For the two- 
bay, multi-story frame shown in Fig. 1lc, one needs to express the rotations 
of the joints in terms of three quantities, say 41, Mz and 63, so that succeed- 
ing equations involve only one additional unknown. The cut is introduced at 
joint 14. The conditions of continuity and equilibrium at the cut may then be 
expressed as 
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+ + = O 
~ = + + 628,50 (28) 
Mi, = + + = O 


where the c’s are again numerical constants. Therefore, the criterion for a 
natural frequency is 


814,13 


C, 


C3, 


C3; 


It is seen that, for this problem, it becomes necessary to evaluate a third 
order determinant. In general, the order of the determinant to be evaluated 
is equal to the number of unknowns used to determine the rotations of the 
joints. 


Simplified Approach 

Even if only two unknowns are needed to express the rotations of the joints, 
the computational work required to calculate the exciting moment and the dis- 
continuity of slope at the cut may get fairly involved. One may remarkably 
simplify this work by determining the desired quantities as a linear combina- 
tion of partial solutions. For illustration, refer again to the frame shown in 
Fig. lla. First, assume that 61 = 1.00 and that @) = 0. If now Eq. (23a) is 
applied to joint 1, it will be found that the rotation of joint 3 will be the only 
unknown in the resulting expression. Determine this unknown and, proceed- 
ing from joint to joint across the frame in the order described previously, 
evaluate the rotations of the remaining joints and designate them by @'. Cal- 
culate also the magnitude of the exciting moment corresponding to the rota- 
tions @' and designate it by M,. Next, assume that 6; = 0 and that 69 = 1.00 
and, in the same manner, calculate the rotations of the joints and the corre- 
sponding exciting moment. Denote these quantities by and M,', respective- 
ly. Then, the actual rotation at a joint j is 


= 66 + (30) 


For the frame shown in Fig. 11a, the total discontinuity of slope at the cut is 


- O65 = (Oss - + Of, = 0, (31) 


and the total exciting moment is 


Mi, = + Mite,. (32) 
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|c,, = 0. (29) 
| 


The coefficients of @; and @2 in these expressions correspond to the coeffi- 
cients c in Eq. (26b). The direct combination of the foregoing partial effects 
is justified by the fact that the differential equation for steady-state forced 
vibration is linear. 


Outline of Procedure 

The details of the main procedure can now be summarized as follows. 
Some of the steps already described are repeated here for the sake of com- 
pleteness. 


1) Assume that the frame is cut at a suitable joint. Then for some mem- 
ber of the frame, say member r, assume a value of A, (See Eq. 5). This is 
equivalent to assuming a frequency of vibration w. 


2) From Eq. (17) compute the values of A for the remaining members of 
the frame, and from the table in the Appendix, obtain the appropriate values 
of K and kk for each member. 


3) The values of K and kk, as given in the Appendix, are expressed in 
terms of Ej1j/Lj, where the subscript j refers to the particular member con- 
sidered. Next express these quantities in terms of the EI/L of the reference 
member r, by multiplying them by the dimensionless factor a; (See Eq. 18). 


4) Compute the sum of the stiffnesses of the members meeting at the var- 
ious joints and record them on a diagram of the frame. On the same diagram, 
record also the product of the stiffness and the carry-over factor for each 
member. These quantities must be expressed in terms of the EI/L of the 
same reference member r. A convenient scheme for arranging the computa- 
tions is shown in the illustrative example presented in the next section. 


5) Choose the unknowns in terms of which the distortions of the frame will 
be expressed. In general, the number of unknowns that must be selected is 
equal to the number of the main longitudinal members in the frame. Consider 
the first of these unknowns equal to unity and the other equal to zero. By 
working across the structure, in the manner described previously, compute 
the rotations of the joints. Denote these by 6’. Compute also the exciting 
moment at the cut and denote it by M,. 


6) Repeat step 5, taking the second unknown equal to unity and the others 
equal to zero. Denote the resulting rotations by 6" and the corresponding 
exciting moment by Mz. In general, steps 5 and 6 must be repeated as many 
times as there are unknowns in the expressions for the distortions of the 
system. 


7) Determine the total discontinuity of slope and the total exciting moment 
at the cut, and evaluate the determinant of the coefficients of the unknowns. 


8) Repeat steps 1 through 7 for a number of values of Ar and plot the vari- 
ation of the determinant evaluated in step 7 as a function of A;. The values 
of Ay for which the determinant becomes equal to zero correspond to the de- 
sired natural frequencies of the frame. 


Illustrative Example 

As an illustration of the application of the method of continuous closed 
frames, we shall calculate the first eight natural frequencies of symmetrical 
vibration of the symmetrical frame shown in Fig. 12. All columns are con- 
sidered to be hinged at the base. 

Since the frame is symmetrical, for symmetrical natural modes of vibra- 
tion the joints of the frame, even though free from external restraining 
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forces, will not translate during free vibration. Accordingly, the procedure 
described is applicable to this problem. In the following, only one half of the 
frame will be considered, The effect of symmetry will be taken into account 
by using for the girders of the central bay modified stiffnesses K® instead of 
the usual stiffnesses K. 

The dimensions of the frame and some additional data pertinent in the an- 
alysis are assembled in Table 2. Member (2) of the frame is taken as the 


1 


ratio of A 
the sake 
computations are given for a value of ho = 2.20, which corresponds toa cir- 


2 
cular frequency of vibration w= 2.20 
L 


The values of a for the various members are recorded in parentheses on 
the diagram in Fig. 13a. The numbers above these values give the stiffnesses 
of the various members, while those below the values of a give the products 
of the stiffnesses and the carry-over factors. These quantities are expressed 
in terms of the EI/L of the member to which they refer, and they are obtained 
directly from the Appendix. 

On the diagram in Fig. 13b, the number without parentheses opposite each 
joint represents the total stiffness of the members framing into that joint. 
These stiffnesses are expressed in terms of the EI/L of the reference mem- 
ber (2). Thus, the total stiffness of the members connecting at joint 4 is 


Ky = + 0.9188(1.167) + 3.943(0.8333) + | = 


r 
0.9294 


[tes | | ons [eres | 


reference member. Columns (5) and (6) of Table 2 give, respectively, the 
y/h2 2 (Eq. 17) and the values of aj (Eq. 18) for each member. For 


TABLE 2 


0.6700 | 0.6700 | 2.08 2.08 


0.50 


brevity, only one cycle of the procedure is presented. Detailed 


735-26 


| 


€ YOs JO 


@) 


02°2 


‘WAS 
> 
(eee!) 
| 


Leee=,"w 
1/13 6692 =,"W 


= 
“ 68 
8B 
<2 
—8 
| 
evo'z 
(€€€9"0) (00°!) ————¢ ro) 
S 
= 
age 252 
Lv0'2 
09" (19:0) 
135-27 


€ 4O NOILONNS 
v sv JO NOILVINWA 


m 
+ 
m 
= 
> 
z 


(29 


= 
“ 
54" 
a 
@ 
> 
f 5 
135-28 


The number at the middle of each member represents the product of the stiff- 
ness and the carry-over factor for that member; these values are also ex- 
pressed in terms of the EI/L of the reference member (2). 

The rotations of the joints are expressed in terms of 9; and @9. The 
frame is assumed to be cut at joint 5 (and joint 5'). First, consider that 
61 = 1.00 and that @2 = 0. Applying Eq. (22) to joints 1 and 2, one finds that 


These values are shown in parentheses in Fig. 13b, and they are recorded on 
the diagram as they are computed. Applying Eq. (23a) successively to joints 
3, 4, and 6, one obtains 


6.001(-1. 799) +1. 426(1. 00) + 52333(0) 12630, 


1.226 


+ 5 +533(-1.799) +0 5.637, 


1.702 


5,6 3.5 8.422, 


The discontinuity in slope at joint 5 corresponding to these distortions is 
' ' 
5,3 95 6 = 16.052 


and the corresponding exciting moment is 


i, = + 1.228(-1.799) + 1.052(-8.422) + 3.555.631) 


EI 
M; = 20.99 L 


Next, consider that 9; = 0 and that 69 = 1.00. The rotations are computed 
in the same manner and the results are recorded on the diagram below the 
values for 6’. 

The discontinuity in slope and the exciting moment at joint 5 correspond- 
ing to these distortions are 


=z 18.2k2 
EI 
= 35.97 


The total discontinuity in slope and the total exciting moment are, 
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2.566(1.00 
3.768(0 


respectively, 


= 16.052 + 18.2he 6, 


5,3 %5,6 


Er = 20.99 + 53.97 . 


The determinant of the coefficients of 61 and @9 is 


16.052 18.242 
4.(6, = = 52.9. 
26.99 33-97 


Since this is different from zero, the assumed value of Az = 2.20 does not 
correspond to a natural frequency. By repeating this procedure for several 
values of Ag and by evaluating, in each case, the resulting determinant, the 
curve given in Fig. 14 was obtained. The values of A» corresponding to the 
natural frequencies of the frame are recorded on the figure. The fundamen- 
tal natural frequency, in radians per second, is 


(2.163)* 
n 


It can be shown that, for the range of frequencies considered, there are no 
natural frequencies corresponding to modes of partial vibration. 


CONCLUDING REMARKS 


A numerical method has been presented for calculating the undamped natu- 
ral frequencies of flexural vibration of continuous beams on rigid supports 
and of rigid-jointed, plane frameworks for which the joints do not displace. 
The members composing the structures have been considered to be of uni- 
form cross section and mass per unit of length. 

In this discussion, the effect of permanent axial forces has been neglected. 
This effect which may be important when the magnitude of the external load 
on the structure is a sizeable fraction of the critical buckling load, may be 
taken into account by using modified stiffness and carry-over factors which 
include the influence of the axial forces. These quantities may conveniently 
be expressed in terms of the dimensionless parameter A and the ratio P/Po, 
where P is the axial compressive or tensile force and Pp is the fundamental 
buckling load of the member considered, assuming its ends to be simply sup- 
ported. For the limiting case of A = 0, the foregoing values of stiffness and 
carry-over factor reduce to those for an axially loaded bar which does not 
vibrate. Thus, the problem of investigating the elastic instability of continu- 
ous beams and frames becomes a special case of the more general problem 
of calculating the natural frequencies of flexural vibration of structures in 
which the component members are subjected to axial forces. (15) 

The procedure may be used also to calculate natural frequencies of beams 
having variable cross section or variable mass per unit of length. A further 
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extension of the method can be made to the case of continuous beams on an 
elastic foundation of the Winkler type. To accomplish this, it is only neces- 
sary to tabulate the pertinent quantities of stiffness and of the product of the 
stiffness and the carry-over factor. 

An obvious extension of the method is in the analysis of the steady-state 
forced vibration of continuous beams and frames subjected to harmonically 
varying forces. For a discussion of this problem the reader is referred to 
reference (15). 

The extension of the method to the case of continuous beams on flexible 
ey and of simple frames with sidesway has been described else- 
where.(16) The method has been applied also to the case of rectangular 
plates which are simply supported along two opposite edges and, in one direc- 
tion, are continuous over rigid supports transverse to the simply supported 


edges, (15 
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APPENDIX 


TABLE I 
DIMENSIONLESS COEFFICIENTS OF K AND kK 
FOR A UNIFORM BAR FIXED AT THE FAR END 


The coefficients are given as a function of the dimensionless 
parameter 4, defined by Eq. (5) in the text. 


BSAAT FINES FURNES SESES 


ke 
l. 
Le 
1. 
l. 
1 
1. 
1. 
l. 
ae 
p 
1 
x 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
Le 
1 
1. 
Ze 
1. 
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L L 4 L 

: (6) 4.0000 | 2.0000 1.40 | 43.9632 | 2.0277 1.80 | 4.8982 | 2.0768 
0.20 | 4.0000 | 2.0000 1.41 | 3.9621 | 2.0285 3.8959 | 2.0785 
0.40 | 3.9998 | 2.0002 1.42 | 43.9610 | 2.0293 3.8935 | 2.0803 
0.60 | 4.9988 | 2.0009 1.43 | 3.9599 | 2.0301 3.8911 | 2.0821 
0.70 | 4.9977 | 2.0017 1.44 | 3.9587 | 2.0310 3.8886 | 2.0840 
0.75 | 3.9970 | 2.0023 3.9576 | 2.0319 3.8862 | 2.0859 
0.80 | 3.9961 | 2.0029 3.9564 | 2.0328 3.8836 | 2.0878 
0.85 | 3.9950 | 2.0037 3.9552 | 2.0337 3.8810 | 2.0898 
0.90 | 3.9937 | 2.0047 3.9539 | 2.0346 4.8784 | 2.0917 
0.95 | 3.9922 | 2.0058 3.9527 | 2.0356 3.8758 | 2.0938 
1.00 | 4.9905 | 2.0072 3.9514 | 2.0366 4.8731 | 2.0958 
1.02 | 4.9897 | 2.0077 3.9500 | 2.0376 4.8703 | 2.0979 ! 
1.04 | 3.9888 | 2.0084 3.9487 | 2.0366 3.8675 | 2.1000 
1.06 | 4.9879 | 2.0090 3.9473 | 2.0396 3.8646 | 2.1022 
1.08 | 4.9870 | 2.0097 3.9459 | 2.0407 3.8617 | 2.1044 
1.10 | 4.9860 | 2.0105 3.9445 | 2.0418 3.8588 | 2.1067 
1.12 | 4.9850 | 2.0113 3.9430 | 2.0429 4.8558 | 2.1089 
1.14 | 3.9839 | 2.0121 4.9415 | 2.0440 3.8528 | 2.1113 
1.16 | 4.9827 | 2.0130 3.9400 | 2.0451 4.84907 | 2.1136 
1.18 | 4.9815 | 2.0139 3.9385 | 2.0463 3.8465 | 2.1160 - 
1.20 | 3.9802 | 2.0149 4.9369 | 2.0475 3.8433 | 2.1184 e 

.|1.21 | 3.9795 | 2.0154 3.93535 | 2.0487 4.8401 | 2.1209 
1.22 | 4.9788 | 2.0159 4.9436 | 2.0500 3.8368 | 2.1234 
1.23 | 3.9781 | 2.016% 4.9319 | 2.0512 3.8334 | 2.1260 
1.24 | 3.9774 | 2.0170 3.9302 | 2.0525 3.8300 | 2.1286 
1.25 | 3.9767 | 2.0175 3.9285 | 2.0536 4.8265 | 2.1312 
1.26 | 4.9759 | 2.0181 4.9267 | 2.0552 4.8230 | 2.1339 
1.27 | 4.9751 | 2.0187 4.9249 | 2.0565 3.8195 | 2.1366 
1.28 | 3.9743 | 2.0193 3.9231 | 2.0579 3.8158 | 2.1394 
1.29 | 3.9735 | 2.0199 3.9212 | 2.0593 4.8121 | 2.1422 
1.30 | 3.9727 | 2.0205 4.9193 | 2.0608 3.8084 | 2.1451 
1.31 | 3.9718 | 2.0212 3.9174 | 2.0622 3.8046 | 2.1480 
1.32 | 3.9709 | 2.0218 3.9154 | 2.0637 3.8007 | 2.1510 
1.33 | 3.9700 | 2.0225 3.9134 | 2.0653 4.7968 | 2.1540 
1.34 | 3.9691 | 2.0232 3.9113 | 2.0668 3.7928 | 2.1570 
1.35 | 3.9682 | 2.0239 3.9092 | 2.0684 3.7887 | 2.1601 
1.36 | 4.9672 | 2.02h6 3.9071 | 2.0700 3.7846 | 2.1643 
1.37 | 3.9662 | 2.0254 3.9049 | 2.0717 4.7805 | 2.1665 
1.38 | 4.9652 | 2.0261 4.9027 | 2.0743 4.7162 | 2.1697 
1.39 | 3.9642 | 2.0269 3.9005 | 2.0750 3.7119 | 2.1730 


TABLE I (Cont'd) 


2- 
2. 
2 
2 
2 
2 
2 
2 
2 
2. 
2 
2 
2. 
2 
2 
2 
2 
2 
2 
2 
2 
2 
2. 
2. 
2s 
2. 
2. 
2. 
2. 
2. 
2. 
2. 
2 
2 
2 
3 
3 
3 
5 
3 
3 
3 
3 
5 
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. . . . . . . . . 


BR 


4 
L a L L 
K KK == K = kK == L 
2.20 | 4.7675 | 2.176% 3.4865 | 2.3947 3.10 | 2.9509 | 2.82h3 ; 
2.21 | 4.7631 | 2.1798 3.4780 | 2.4014 3.11 | 2.9346 | 2.8377 
2.22 | 4.7586 | 2.1832 3.4693 | 2.4082 3.12 | 2.9180 | 2.8513 
2.23 | 3.7540 | 2.1868 3.4605 | 2.4151 3.13 | 2.9011 | 2.8651 
2.24 | 3.749% | 2.1903 3.4516 | 2.4222 3.14 | 2.8841 | 2.879]. 
2.25 | 3.7447 | 2.1939 3.4426 | 2.4293 3.15 | 2.8667 | 2.8933 
2.26 | 4.7399 | 2.1976 3.43354 | 2.4365 3.16 | 2.8491 | 2.9078 
2.27 | 3.7350 | 2.201% 3.4241 | 2.4439 3.17 | 2.8312 | 2.9225 
Bs 2.28 | 4.7301 | 2.2052 3.4147 | 2.4513 3.18 | 2.8131 | 2.937% : 
2.29 | 3.7251 | 2.2090 3.4051 | 2.4589 3.19 | 2.7947 | 2.9526 
2.30 | 4.7200 | 2.2129 4.3954 | 2.4666 3.20 | 2.7760 | 2.9681 : 
2.31 | 3.7149 | 2.2169 3.3855 | 2.474% |] 3.21 | 2.7569 | 2.9836 
2.32 | 4.7096 | 2.2209 3.3755 | 2.4823 3.22 |. 2.7377 | 2.9997 
2.33 | 3.7043 | 2.2250 3.3654 | 2.4904 3.23 | 2.7181 | 3.0160 
2.34 | 4.6989 | 2.2292 3.3551 | 2.4985 3.24 | 2.6982 | 5.0325 
2.35 | 3.6935 | 2.334 3.3447 | 2.5068 3.25 | 2.6779 | 3.0493 
2.36 | 3.6879 | 2.2577 3.3341 | 2.5152 3.26 | 2.6574 | 3.0663 
2.37 | 3.6824 | 2.2321 3.3234 | 2.5237 3.27 | 2.6366 | 34.0837 
2.38 | 3.6766 | 2.2465 3.3125 | 2.5324 3.28 | 2.6154 | 3.1013 
2.39 | 3.6708 | 2.2510 3.3015 | 2.5412 3.29 | 2.5939 | 3.1193 : 
2.40 | 3.6649 | 2.2555 4.2903 | 2.5501 3.30 | 2.5720 | 3.1375 : 
2.41 | 4.6589 | 2.2601 4.2790 | 2.5592 3.31 | 2.5498 | 3.1561 
2.42 | 4.6528 | 2.2648 3.2674 | 2.5684 3.32 | 2.5272 | 3.1750 
2.43 | 4.6467 | 2.2696 3.2558 | 2.5778 3.33 | 2.5043 | 3.1942 
2.44 | 3.6405 | 2.2744 3.2439 | 2.5872 3.44 | 2.4810 | 3.2138 
2.45 | 3.6341 | 2.2793 3.2319 | 2.5969 3.35 | 2.4573 | 3.2337 
2.46 | 3.6277 | 2.2843 4.2197 | 2.6067 3.36 | 2.4332 | 3.2539 
2.47 | 4.6212 | 2.2894 3.2073 | 2.6166 3.37 | 2.4087 | 3.2746 
2.48 | 4.6146 | 2.2945 3.1948 | 2.6267 3.38 | 2.3839 | 4.2955 
2.49 | 4.6079 | 2.2997 4.1820 | 2.6369 3.39 | 2.3586 | 4.3169 
2.50 | 3.6011 | 2.3050 3.1691 | 2.6473 3.40 | 2.3329 4386 
2.51 5942 | 2.4104 4.1560 | 2.6579 3.41 | 2.3067 3607 
2.52 5872 | 2.3159 3.1427 | 2.6686 3.42 | 2.2802 3833 
2.53 5801 | 2.3214 4.1292 | 2.6795 5.44 | 2.2531 4062 
2.54 5729 | 2.3270 4.1155 | 2.6906 3.44 | 2.2256 4296 
2.55 5656 | 2.3327 4.1016 | 2.7018 3.45 | 2.1977 4534 
2.56 5582 | 2.3385 3.0875 | 2.7132 3.46 | 2.1692 4776 
2.57 5507 | 2.3444 3.0732 | 2.7248 3.47 | 2.1403 5023 
2.58 5430 | 2.3504 4.0587 | 2.7366 3.48 | 2.1109 5274 
2.59 5353 | 2.3564 3.0440 | 2.7485 3.49 | 2.0809 5541 
2.60 5275 | 2.3625 3.0290 | 2.7607 3.50 | 2.0505 5792 
2.61 5195 | 2.3688 3.0139 | 2.7730 3.51 | 2.0194 6058 
2.62 5114 | 2.3751 2.9985 | 2.7855 3.52 | 1.9879 6329 
2.63 5033 | 2.3815 2.9828 | 2.7983 3.53 | 1.9557 5 
2.64 4O49 | 2.3881 2.9670 | 2.8112 3.54 | 1.9230 | 3 7 
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TABLE I (Cont'd) 


On ow Fu 


FEF EE FREE EE FREE EE FE 
SES SESES FUDES SSQR 


L L L 
55 | 1.8897 3.7174 114.05 | -1.0838 | 6.3997 4.55 | -20.539 25.402 
56 | 1.8558 3.7467|1 4.06 | -1.1891 | 6.4978 4.56 | -22.088 26.941 
57 | 1.8213 4.7766 4.07 | -1-2977 | 6.5991 4.57 | -23.841 28.672 
58 | 1.7861 4.8071 || 4.08 | -1.4096 | 6.7036 4.58 | -25.805 50.635 
59 | 1.7503 3.83811] 4.09 | -1.5249 | 6.8115 4.59 | -28.060 42.879 
60 | 1.7138 4.8698|| 4.10 | -1.6439 | 6.9230 -30.662 35.469 
61 | 1.6766 3.9022|] 4.11 | -1.7668 | 7.0383 -34.697 38 .4ug2 
62 | 1.6387 3.9352|] 4.12 | -1.8937 | 7.1575 -37.282 42.066 
63 | 1.6001 4.9689 114.13 | -2.0247 | 7.2609 -41 46.355 
64 | 1.5607 4.0033]1] 4.14 | -2.1603 | 7.4086 -46.849 51.599 
65 | 1.5206 4 .0364}14.15 | -2.3004 | 7.5410 -53.407 58.155 
66 | 1.4797 4.07K2 114.16 | -2.4455 | 7.6781 -61.849 66.585 
67 | 1.4380 4.1108 }14.17 | -2.5957 | 7.8203 -73.102 17.826 
68 | 1.3954 4.1482114.18 | -2.7513 | 7.9679 -88 .850 93.561 
69 | 1.3520 4.1864]}14.19 | -2.9127 | 8.1212 -112.46 117.16 
70 | 1.3078 | 4.2254]14%.20 | -3.0802 | 8.2805 -151.79 156.48 
71 | 1.2626 4.2653]] 4.21 | -5.2540 | 8.4461 -230.36 235.04 
72 | 1.2165 4.3060} ]4.22 | -3.4347 | 8.6183 -465 470.09 
73 | 1.169% 4.3477} 14.23 | -3.6225 | 8.7977 oo too 
1.1214 4.39021] 4.2% | -3.8179 | 8.9847 480.58 |-475.05 
75 | 1.0724 4.43381 | | -4.0214 | 9.1797 242.62 |-238.00 
76 | 1.022% 4.4784 26 | -4.2335 | 9.3832 164.51 |-158.91 
77 | 0.97110 | 4.5239 27 | -4.4548 | 9.5958 124.00 |-119.40 
78 | 0.91884 | 4.5705 28 | -4.6859 | 9.8180 100.40 -95.717 
79 | 0.86543 | 4.6182 29 | -4.9274 [10.051 84.501 | -79.934 
80 | 0.81084 | 4.6670 30 | -5.1801 |10.29% 73.219 | -68.666 . 
81 | 0.75504 | 4.7170 31 | -5.4447 110.550 64.757 | -60.218 
82 | 0.69799 | 4.7682 32 | -5.7221 {10.818 58.176 | -53.650 
84 | 0.63965 | 4.8206 33 | -6.0133 |11.100 52.909 | -48.398 
84 | 0.57997 | 4.8744 34 | -6.3194 |11.397 48.600 | -44.104 7 
85 | 0.51891 | 4.929% 35 | -6.6415 [11.710 45.008 | -40.525 
86 | 0.45642 | 4.9858 36 | -6.9809 |12.040 41.967 | -37.499 
87 | 0.39246 | 5.0436 37 | -7.3390 |12.389 | 39.360 | -34.907 
88 | 0.452697 | 5.1029 38 | -7.71L74% |12.757 37.100 | -32.661 
89 | 0.25990 | 5.16386 49 | -8.1180 {13.148 45.122 | -30.698 = 
90 | 0.19119 | 5.2262 ko | -8.5427 113.564 43.375 | -28.966 
91 | 0.12079 | 5.2902 4) | -8.9938 |14.005 41.822 | -27.428 
92 | 0.048644) 5.3560 42 | -9.4739 114.475 30.431 | -26.053 
93 |-0.025344| 5.4235 43 | -9.9858 |14.977 29.179 | -24.816 
94 |-0.10121 | 5.4928 4h 1-10.533 {15.514 28.045 | -23.698 
-0.17904 | 5.5641 45 |-11.119 |16.089 27.014 | -22.682 
-0.25890 | 5.6373 46 |-11.748 [16.708 26.072 | -21.756 
-0.34087 | 5.7126 47 |-12.425 117.375 25.207 | -20.907 
-0.42504 | 5.7900 48 |-13.156 {18.096 24.411 | -20.127 
-0.51150 | 5.8696 4g |-13.948 |18.877 23.676 | -19.408 
-0.60034 | 5.9516 50 |-14.808 [19.727 22.994 | -18.7%3 
-0.69164 | 6.0460 51 |-15.747 |20.655 22.361 | -18.126 
-0.78553 | 6.1229 52 |-16.77% 21.771 | -17-553 
-0.88210 | 6.212% 53 |-17.904 {22.790 21.219 | -17.018 
-0.98148 | 6.3046 54 |-19.152 |24.027 20.702 | -16.518 
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20.217 
19.761 
19.3531 
18.925 
18.542 


18.178 
17.833 
17.505 
17.193 
16.895 


16.612 
16.541 
16.061 
15-833 
15.595 


15.367 
15 
14.937 
14.734 
14.538 


14.350 
14.166 
13.992 
15.623 
13.658 


13.500 
13.346 
13.197 
13.052 
12.911 


12.643 
12.514 
12.388 
12.266 


12.147 
12.031 
11.918 
11.808 
11.701 


11.595 
11.493 
11.392 
11.294 
11.198 


11.104 
11.011 
10.921 
10.832 
10.745 


. 


. . . . . 
FUBES SELES SE 


@ 


= 


5. 
5. 
5. 
5. 
5. 
5. 
5. 
5. 
5. 
5. 
5. 
5. 
5. 
5. 
5. 
5. 
5. 
5. 
5. 
5. 
5. 
5. 
5. 
5. 
5. 
5. 
5. 
5. 
5. 
5 
3 
5 
5 
5 
5. 
5 
5 
5. 
5. 
5. 


“fs 
\ 
; L L L L L L 

05 -16.050 5.55 |10.660 -7-5873 || 6. 7.5233 | -6.2553 
06 -15.611 5.56 |10.576 -7.5312 || 6. 7.4711 | -6.2504 
OT -15.199 5.57 10.494 -7.4768 || 6. 7.4190 | -6.2461 
08 -14.811 5.58 110.413 -7.4242 || 6. 7.3670 | -5.2625 
“14.445 5.59 110.334 -7.3732 || 6. 7.3151 | -6.2396 
10 -14.099 5.60 110.256 -7.3239 || 6. 7.2632 | -6.2573 
11 -13.772 5.61 |10.179 -7.2761 || 6. 7.2113 | -6.2356 
12 -13.462 5.62 |10.103 -7.2298 || 6. 17-1595 | -6.2346 
13 -13.169 5.63 |10.029 -7.1850 || 6. 7.1077 | -6.2342 
14 -12.890 5.64 | 9.9555 | -7-1416/| 6. 7.0559 | -6.2344 

15 -12.625 5.65 | 9.86832 | -7.0996 || 6. 7-O041 | -6.2353 
16 -12.373 9.8120 -7.0590 6. 6.9522 -6.2368 
17 -12.133 9.7418 -7-0196 6. 6.9003 -6.2389 
18 -11.904 9.6726 | -6.9815 || 6. 6.8484 | -6.2416 
19 -11.685 9.6042 | -6.9447 1) 6. 6.7964 | -6.2450 
20 -11.477 9.5368 | -6.9091 || 6. 6.7444 | -6.2490 
21 -11.277 9.4702 | -6.8747 || 6. 6.6922 | -6.2536 
22 -11.086 9.4044 | -6.8414 |] 6. 6.6400 | -6.2589 
23 -10.903 9.339% | -6.8092 || 6. 6.5876 | -6.2648 
24 -10.728 9.2752 | -6.7781 || 6. 6.5351 | -6.2713 
25 -10.560 9.2116 | -6.7481 || 6. 6.4825 | -6.2785 
26 -10.399 9.1488 | -6.7192 || 6. 6.4297 | -6.2862 
27 -10.245 9 .0866 -6.6912 6. 6.3768 -6.2946 
28 -10.096 9.0251 | -6.6643 || 6. 6.3237 | -6.3037 
29 - 9.9534 8.9642 | -6.6364 || 6. 6.2704 | -6.3134 
40 - 9.8160 8.9038 | -6.6134 || 6.30 | 6.2169 | -6.3237 
31 - 9.6839 8.8441 | -6.5893 || 6.31 | 6.1631 | -6.3347 
42 - 9.5567 8.7848 | -6.5662 || 6.32 | 6.1092 | -6.3463 
43 - 9.4542 8.7261 | -6.5439 || 6.33 | 6.0550 | -6.3585 
a4 - 9.3162 8.6678 | -6.5226 || 6.34 | 6.0006 | -6.3715 
35 | 12.775 | - 9.202% 8.6101 | -6.5021 || 6.35 | 5.9459 | -6.3650 
46 - 9.0926 8.5527 | -6.462h || 6.36 | 5.8909 | -6.3993 
47 = 8.9868 8.4958 -6.4636 6.37 5 8356 -6.4142 
- 8.8846 8.4393 | -6.4457 || 6.38 | 5.7800 | -6.4298 
- 8.7860 8.3832 | -6.4285 || 6.39 | 5.7241 | -6.4461 
- 8.6908 Bo | © | 5.6679 | -6.4630 
a 8.5988 1 5.6113 -6.4807 
- 8.5099 2 | 5.554% | -6.4990 
> 8.4240 b) 5.4970 -6.5181 
- 8.3409 4b | 5.4393 | -6.5379 
- 8.2606 5 | 5.3612 | -6.5584 
- 8.1829 6 | 5.3226 | -6.5796 
8.1078 7 5.2636 -6.6016 
@ 8.0351 7.8925 -6.3086 6.48 5.2042 -6.6243 
- 71-9647 7.8392 | -6.2989 || | 5.1442 | -6.6478 
- 17-6966 7.7862 | -6.2899 || | 5.0838 | -6.6720 
- 7.8307 7.7333 | -6.2617|| 6.81 | 5.0229 | -6.6971 
- 7.7668 7.6806 | -6.27%1 || 6.fI2 | 4.9615 | -6.7229 
= 7.7051 7.6280 -6.2671 6. 4 4.8995 -6.7496 
- 7.6452 7.5756 | -6.2609|| 6.84 | 4.6369 | -6.7770 

735-36 


TABLE I (Cont'd) 


wy 


OON AM FUNrFO O 


BES FERES 


FEE E 


-64 
.65 
-66 
-67 
-68 
69 
70 
T2 
73 
-75 
-76 
-78 
79 
80 
81 
82 
83 
84 
85 
86 
-88 
89 
90 
gl 
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+95 
96 
OT 
-98 
99 


6 
6 
6 
6 
6 
6 
6. 
6. 
6. 
6. 
6. 
6 
6 
6 
6 
6 
6. 
6. 
6. 
6. 
6. 
6 2. 
6. 2. 
6 2. 
6 
6. 2. 
6. 2. 
6. Li 
6. 
6. 
6 
6 Zs 
6 
6 1. 
6 
6. 1. 


— 
L L L L 
6.55 | 4.7738 | -6.8053]|7.00 | 0.90203 | -9.2964//7.45]| -10.007] -18.989 
6.56 | 4.7101 | -6.8344]]7.01 | 0.77860 | -9.3919]/7-46| -10.516] -19.471 
6.57 | 4.6458 | -6.8644]1]7.02 | 0.65256 | -9.4901]/7.47| -11.052] -19.980 
6.58 | 4.5808 | -6.8954 ; -11.616 | -20.517 
6.59 | 4.5152 | -6.9270 ‘. -12.209 | -21.084 
6.60 | 4.4489 | -6.9597 J ‘ -12.836 | -21.683 
6.61 | 4.3818 | -6.9933 ‘ -13.498 | -22.319 
6.62 | 4.3141 | -7.0278 -14.200 | -22.993 
6.63 | 4.2457 | -7.0632 R -14.944 | -23.710 
1764 | -7.0997 : -15.734 | -24.473 
1064 | -7.1371 x -16.576 | -25.287 
0356 | -7.1755 -17.474 | -26.158 
9639 | -7.2150 ; -18.434 | -27.091 
8914 | -7.2555 : -19.464 | -28.094 
8179 | -7.2971 -20.571 | -29.173 
7436 | -7.3398 -21.764 | -30.339 
6683 | -7.3836 -23.054 | -31.602 
5920 | -7.4286 24.454 -32.974 
5147 | -7.4747 -25.977 | -34.470 
4364 | -7.5220 -27.643 | -36.108 
3570 | -7.5706 -29.470 | -37.908 
2765 | -7.6204 -31.486 | -39.895 
1949 | -7.6715 -33.720| -42.101 
1121 | -7.7240 -36.210 | -44.564 
0281 | -7.7778||7 -39.003 | -47.329 
9428 | -7.83301)7 -42.159 | -50.457 
8563 -7.8896 “45.754 -54.024 
7684 | -7.9476117 -49.886 | -58.128 
6791 | -8.0072]|7 | -54.687 | -62.900 
5885 | -8.0683|17 -60.333 | -68.518 
4963 | -8.1310117 -67.070 | -75.227 
4027 | -8.1953|17 -75-250 | -83.379 
4075 | -8.2613117 -85.392 | -93.492 
2107 | -8.3290] 17 -98.301 | -106.37 
1123 | -8.3985/|7 -115.29 | -125.33 
0121 | -8.4698} | 7 -138.66 | -146.67 
9101 | -8.5430 -172.84 | -180.82 
8064 | -8.6181 -227.60 | -235.55 
7007 | -8.6952 -329.54% | -337.46 
5930 | -585.85 | -593.75 
4834 | -8.8556 2441.7 2449.6 
3716 | -8.9391 1164.5 1156.7 
2577 | -9.02%9 
1415 | -9.1129 | 
0230 | -9.2034 
735-37 


Term expires October, 1955: 
CHARLES B. MOLINEAUX 
MERCEL J. SHELTON 
A. A. K. BOOTH 
CARL G. PAULSEN 
LLOYD D. KNAPP 
GLENN W. HOLCOMB 
FRANCIS M. DAWSON 


OFFICERS FOR 


Term expires October, 1955: 
ENOCH R. NEEDLES 
MASON G. LOCKWOOD 


WALTER L. HUBER 


EXECUTIVE SECRETARY 
WILLIAM H. WISELY 


ASSISTANT SECRETARY 
E. L. CHANDLER 


DEFOREST A. MATTESON, JR. 
Editor of Technical Publications 


GLENN W. HOLCOMB 
ERNEST W. CARLTON 


AMERICAN SOCIETY OF CIVIL ENGINEERS 


1955 


PRESIDENT 
WILLIAM ROY GLIDDEN 


VICE-PRESIDENTS 


DIRECTORS 


Term expires October, 1956: 
WILLIAM 5S. LaLONDE, JR. 
OLIVER W. HARTWELL 
THOMAS C. SHEDD 
SAMUEL B. MORRIS 
ERNEST W. CARLTON 
RAYMOND F. DAWSON 


'PAST-PRESIDENTS 
Members of the Board 


PROCEEDINGS OF THE SOCIETY 


HAROLD T. LARSEN 


COMMITTEE ON PUBLICATIONS 
SAMUEL B. MORRIS, Chairman 


JEWELL M. GARRELTS, Vice-Chairman 
OLIVER W. HARTWELL 


DON M. CORBETT 


Term expires October, 1956: 
FRANK L. WEAVER 
LOUIS R. HOWSON 


DANIEL V. TERRELL 


CHARLES E. TROUT 


ASSISTANT TREASURER 
CARLTON S. PROCTOR 


Manager of Technical Publications 


Assoc. Editor of Technical Publications 


Term expires October, 1957: 
JEWELL M. GARRELTS 
FREDERICK H. PAULSON 
GEORGE 5S. RICHARDSON 
DON M. CORBETT 
GRAHAM P. WILLOUGHBY 
LAWRENCE A. ELSENER 


TREASURER 


PAUL A. PARISI 


— 


